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Abstract: This review examines the fundamentals and challenges in engineering/understanding the thermostability of bio- 
logical systems over a wide temperature range (from the cryogenic to hyperthermic regimen). Applications of the bio- 
thermo stability engineering to either destroy unwanted or stabilize useful biologicals for the treatment of diseases in mod- 
ern medicine are first introduced. Studies on the biological responses to cryogenic and hyperthermic temperatures for the 
various applications are reviewed to understand the mechanism of thermal (both cryo and hyperthermic) injury and its 
quantification at the molecular, cellular and tissue/organ levels. Methods for quantifying the thermophysical processes of 
the various applications are then summarized accounting for the effect of blood perfusion, metabolism, water transport 
across cell plasma membrane, and phase transition (both equilibrium and non-equilibrium such as ice formation and glass 
transition) of water. The review concludes with a summary of the status quo and future perspectives in engineering the 
thermostability of biological systems. 
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1. INTRODUCTION 

Efforts on engineering the thermostability of biological 
systems could be traced back to thousands of years ago ac- 
cording to the records of early civilizations. For example, 
high temperatures have been used to treat a variety of dis- 
eases including back pain in the traditional Chinese medicine 
[1-3]; the ancient Hindus heated metal bars to stop bleeding 
[4]; in 500 B.C., the Greek physician Parmenides believed 
that if he could create fever, he could cure all illness [5]; 
dried mummies have been found in Egypt and other places to 
be well-preserved at ambient temperatures for thousands of 
years [6]; and the use of cold was recommended by Hippo- 
crates to reduce bleeding and swelling [7]. Nothing speaks 
for the unlimited potential of temperature in medicine better 
than an aphorism of Hippocrates [8,9]: "Those diseases 
which medicines do not cure, iron (author's note: iron means 
scalpel here) cures; those which iron cannot cure, fire 
(author's note: fire means temperature here) cures; and those 
which fire cannot cure, are to be reckoned wholly incurable." 

Indeed, temperature (either heat or cold) is still one of the 
most important tools for the treatment of diseases in modern 
medicine. Recent research on engineering the (thermo) sta- 
bility of biologicals at various temperatures for biomedical 
applications has been focused on either destroying un- 
wanted/diseased cells/tissues or stabilizing cells/tissues for 
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future use in the treatment of diseases as summarized in Fig. 
(1). The use of abnormally low and high temperatures to 
locally destroy unwanted tissues is called cryosur- 
gery/cryoab lation/cry other apy and thermal surgery /thermal 
ablation/thermal therapy/thermotherapy (including hyper- 
thermia for which the temperature is usually between 39 and 
45 °C), respectively. For thermal stabilization (or biopreser- 
vation) applications, different terminologies including cryo- 
preservation, hypothermic preservation and dry or lyopreser- 
vation are used for temperatures below -80 °C, at ~ 4 °C, and 
at ambient temperatures, respectively. There is a significant 
overlap in terms of temperature range between hypothermic 
and lyopreservation. However, they are quite different since 
hypothermic preservation is designed for preservation of 
fully hydrated tissues/organs in a short period usually from 
days to weeks while lyopreservation is aimed at stabilizing 
biologicals in the long term (i.e., up to years) in a dehydrated 
state. More detailed discussions of the different applications 
are given below. 

1.1. Thermal Destruction 

On the thermal destruction side, whole body hyperther- 
mia (usually with brain cooling to avoid thermal damage to 
the extremely heat susceptible cerebral tissue [10,11]) with 
temperatures usually between 39-45 °C has been used for 
many years as an adjuvant treatment to enhance the effect of 
chemo and radio therapy of malignant diseases (i.e., cancer) 
[12-16]. More recently, thermal therapy and cryosurgery 
have been developed as single treatment modalities for vari- 
ous benign and malignant diseases [17-22]. With the ad- 
vancement of modern surgical techniques such as laparo- 
scopy, both thermal therapy and cryosurgery can be per- 
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Fig. (1). A summary of typical biomedical applications for which the goal is to engineer the thermostability (for either destruction or stabili- 
zation) of biological systems from cryogenic to hyperthermic temperatures. 
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Fig. (2). Sketches of a typical cryoprobe (A) designed based on the Joule-Thomson effect for destroying tissue by freezing and a typical 
thermal probe (B) using energy carried either in an electromagnetic (e.g., radio frequency, microwave, or laser) or acoustic wave field 
(HIFU, high intensity focused ultrasound) to destroy tissue by heating. 



formed in a minimally invasive manner to reduce the patient 
suffering and hospital stay during and post surgical operation 
[23-25]. Consequently, they are becoming increasingly 
popular and have been studied in essentially all surgical sub- 
specialties as a minimally invasive alternative to the conven- 
tional radical surgical intervention [26-28]. 

1.1.1. Thermal Destruction at Cryogenic Temperatures 
(Cryotherapy): Cooling Approaches 

In cryosurgery, unwanted or diseased tissue is destroyed 
by freezing it to subzero temperatures (usually below -20 °C) 
typically using one or multiple cryoprobes (up to a few mil- 
limeters in diameter, Fig. 2A). The tissue is frozen through 
direct contact with the cryoprobe which is cooled to a cryo- 
genic temperature by circulating cryogenic media at a tem- 
perature usually below -100 °C inside the cryoprobe lumen 
[29,30]. Liquid nitrogen (-196 °C) is one of the commonly 
used cryogenic media, particularly in cryoprobes designed 
before the 1990s. Gases such as argon, carbon oxide, and 
nitrous oxide at cryogenic temperatures obtained by throt- 
tling (i.e., from high to low pressure) as a result of the Joule- 
Thomson effect [31] are becoming increasingly popular. 
Liquid nitrogen has been replaced with these gases in more 
recent designs of cryoprobes for better control of the probe 
temperature. An additional advantage of cryoprobes de- 
signed based on the Joule-Thomson effect is that the circu- 
lating gas for freezing in the cryoprobe can be switched to a 



different gas such as helium, whose temperature actually 
increases after throttling, due to a negative Joule-Thomson 
coefficient (i.e., |i JT defined in Fig. 2A) of the gas at around 
room temperature. This coefficient is positive for argon at 
room temperature or below. This feature is utilized to 
achieve fast, active (instead of slow, passive) thawing of the 
diseased tissue after freezing, which is important for accurate 
control of the frozen tissue volume. Other gas mixtures 
(mainly Freon) at cryogenic temperatures (usually higher 
than -130 °C) obtained utilizing a refrigeration cycle system 
essentially like that of a household refrigerator/freezer, are 
also used in some cryoprobe designs. 

1.1.2. Thermal Destruction at Hyperthermic Temperatures 
( Thermoth erapy ) : Heating Approach es 

In thermal therapy, unwanted or diseased tissue is ther- 
mally ablated or destroyed by locally heating the tissue to 
above 50 °C for minutes or even seconds typically using one 
(for removing a small unwanted tissue volume) or an array 
(for removing a large unwanted tissue volume) of thermal 
probes (also up to a few millimeters in diameter, Fig. 2B). 
Unlike the cryoprobes for cryosurgery that must have a cold 
surface to freeze tissue, a thermal probe for thermal therapy 
is not necessary to be hot and even not necessary to be in 
direct contact with the tissue, particularly when laser or high 
intensity focused ultrasound (HIFU) is used as the energy 
source [32-37]. This is because it is the electromagnetic or 
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acoustic wave field applied on the tissue through the thermal 
probe that heats up and thus destroys the tissue. Moreover, a 
cooling mechanism might be employed to protect the ther- 
mal probe from being overheated when the probe does have 
direct contact with the tissue to be thermally destroyed 
[38,39]. This is because the efficacy of the probe for emitting 
electromagnetic or acoustic wave might be significantly 
compromised at the high temperature necessary for destroy- 
ing unwanted tissue. 

The electromagnetic wave over a wide range of fre- 
quency (or wavelength) including radiofrequency (RF), mi- 
crowave, and laser, has been used for heating. The terms of 
radiofrequency (RF) and microwave have often been used 
interchangeably in the hyperthermic literature for electro- 
magnetic waves with frequency over a broad range from 
several hundred kilohertz (kHz) to several hundred gigahertz 
(GHz). Traditionally, RF is referred as electromagnetic wave 
with frequency from 300 kHz to 1GHz and microwave is 
from 300 MHz to 300 GHz. The corresponding wavelength 
of RF and microwave in air is ~ 1000 - 0.3 m and ~ 1000 - 1 
mm, respectively. Considering the speed of light in water 
(the major building component of most soft tissue) is ap- 
proximately half of that in air (3><10 8 m/s), the wavelength of 
RF/microwave in most soft tissue should decrease by half as 
well. Typical RF/microwave frequencies used for heating 
biological tissue in the hyperthermic literature are around 
2.45 GHz, 915 MHz, 40.68 MHz, 27.12 MHz, and 13.56 
MHz, which are the frequencies allocated by the Federal 
Communications Commission (FCC) for medical use in 
USA. The two highest frequencies (2.45 GHz and 915 MHz) 
are usually (although not always) referred as microwave 
while the rest are called RF in the thermal therapy literature. 

Both conductive (more important at low frequencies) and 
dielectric (more important at high frequencies) heating 
mechanisms contribute to the absorption of electromagnetic 
energy by biological tissue [40-42]. The former is due to the 
electrical resistance of biological tissue and is the major heat 
generation mechanism at low frequencies (e.g., < 1 MHz). 
The latter is due to friction between molecules in tissue as a 
result of the alternating movement of polar molecules (par- 
ticularly the water dipoles in a biological tissue) trying to 
align themselves with the alternating RF/microwave electric 
field. Therefore, dielectric heating is frequency dependent 
and is the major heating mechanism at high frequencies (e.g., 
> 100 MHz), particularly in the nonconductive biological 
tissue. A more comprehensive review of the biological effect 
of RF/microwave can be found elsewhere [40]. 

Lasers (acronym for light amplification by stimulated 
emission of radiation) of infrared, visible light, and ultravio- 
let with a wave length much less than 1 mm (typically -150 
nm - 100 urn in the infrared domain) have also been used for 
thermal destruction. Correspondingly, the wave frequency is 
much higher than 300 GHz but usually less than 10 16 Hz. In 
this case, the tissue is heated by absorbing photons carried in 
the laser wave field [32]. Due to the small wave length, laser 
wave could be scattered significantly along its propagation 
path in biological tissue which is generally inhomogeneous. 
Light scattering can significantly enhance tissue absorption 
of the laser energy on one hand and limit the penetration 
depth of the laser wave in tissue on the other. Generally, the 
heating depth of a single interstitial laser probe is only ap- 
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proximately one tenth of that of an interstitial RF/microwave 
probe under comparable working conditions. Therefore, mul- 
tiple probes or multiple laser shots at different locations in 
the targeted tissue are necessary to destroy a large tissue vol- 
ume. It is interesting and important to note that lasers with 
wavelength from 800 to 1,200 nm can be absorbed much 
more strongly by metals than water [33]. Therefore, lasers in 
this wavelength range can be utilized to specifically destroy 
unwanted tissue by heating metallic (e.g. gold) nanoparticles 
delivered into the tissue [43,44]. 

The acoustic (or mechanical) wave used for heating and 
destroying unwanted tissue is generally referred as high in- 
tensity focused ultrasound (HIFU). The frequency of medical 
acoustic waves (a mechanical wave) is typically from 1 to 3 
MHz. As a result, the wavelength of the HIFU mechanical 
wave is typically less than -1.5 mm since the speed of 
sound in water is ~ 1500 m/s between 0-100 °C. Therefore, 
tissue scattering of the mechanical wave is significant, which 
can enhance its absorption by tissue and at the same time 
limit its penetration depth in tissue. The latter is due to unde- 
sired heating of the surface tissue such as the dermis. Conse- 
quently, multiple HIFU shots at different locations in the 
targeted tissue are necessary to destroy a large diseased tis- 
sue volume. At a low intensity, tissue absorption of acoustic 
energy can be attributed to interstitial mechanical friction as 
a result of the alternating movement in tissue driven by the 
alternating acoustic (mechanical) wave while at a very high 
intensity, cavitation will occur as a result of microbubble 
formation/collapse and its interaction with the acoustic wave 
[35]. A significant enhancement of heat generation has been 
observed when cavitation occurs [36,37]. In addition, the 
generation of shock wave associated with microbubble col- 
lapse during cavitation can destroy tissue mechanically [35]. 

1.2. Thermal Stabilization 

The goal of thermal stabilization (or biopreservation) is 
to preserve important biologicals such as biomolecules, cells, 
tissues, and organs for future use (e.g., transplantation) in 
modern medicine. For short term stabilization, biologicals 
are usually cooled to around ice temperature (hypothermic 
preservation) without ice formation or freezing. The purpose 
is to slow down metabolism and degradation of the biologi- 
cals to extend their shelf life to days and up to weeks after 
explantation [45-47]. Studies on short term hypothermic 
preservation are mainly focused on organs, for which long- 
term stabilization has not been practically successful. For 
long-term stabilization, however, biologicals (e.g., proteins, 
cells and small tissue slices) are usually processed and stabi- 
lized in an amorphous or glassy phase, which can be best 
demonstrated in an extended phase diagram shown in Fig. 
(3). The glassy phase is a thermodynamic ally metastable 
state with an extremely high viscosity and low molecular 
mobility and activity. Therefore, any processes that require 
the diffusion of participating molecules including both de- 
gradative and life supporting processes in a biological sys- 
tem are essentially arrested/suspended in a desired glassy 
phase. Practically, five different approaches including the 
conventional slow-freezing, conventional vitrification using 
a high concentration of cryoprotectant (or cryoprotective 
agent, or CPA for short), low- CPA vitrification, freeze- 
drying, and evaporative/convective drying (or desiccation) 
have been studied to bring biologicals from an aqueous liq- 
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Fig. (3). An illustration on the extended phase diagram of the various approaches to achieve thermal stabilization on biologicals at either 
cryogenic or ambient temperature from an initial liquid state (A) to a final glassy phase (e.g., I, II, III, or IV): The phase diagram is divided 
into four thermodynamic regimens by the liquidus, extended liquidus, solidus, and the glass transition curve; the four regimens are the liquid, 
subcooled liquid, supersaturated liquid, and the glassy phase; cells (in a liquid solution initially) must enter the glassy phase for long-term 
preservation; CPA represents cryoprotectant, and 0 g and 0 m represent the glass transition and melting temperature, respectively. Of note, the 
diagram is not to scale (for example, the melting temperature of pure CPAs is usually below 20 °C). Figure reprinted from reference [53] 
with permission from Xiaoming He (the author). 



uid state at room temperature to a glassy phase at either a 
cryogenic or ambient temperature as shown in Fig. (3). The 
former is usually called cryopreservation while the latter is 
usually called dry or lyopreservation. A more detailed dis- 
cussion of each of the approaches for long-term cell/tissue 
preservation is given in the following sub-sections. 

1.2.1. Thermal Stabilization at Cryogenic Temperature 
(Cryopreservation ) 

Thermal stabilization of biologicals at a cryogenic tem- 
perature can be achieved with two major approaches: the 
conventional slow-freezing with inevitable ice formation and 
the ice-free vitrification [48-53]. For the conventional slow- 
freezing approach, the following steps are typically adopted 
(A— >B— >C— >D— >E— >I in Fig. 3): (1), biologicals such as 
cells in an aqueous solution (state A) are first loaded with 
CPA(s) such as DMSO (dimethylsulf oxide), glycerol, pro- 
pylene glycol (or 1,2-propanediol), and ethylene glycol at a 
concentration usually up to 15 wt% (or up to 2 M, state B); 

(2) , the samples are sub-cooled usually to between -2 and -7 
°C (phase C: note, phase rather than state is used to indicate 
that it is not in equilibrium) to seed ice in the extracellular 
space by touching the sample contained usually in a cryovial 
with a deeply cooled (e.g., in liquid nitrogen) object (C— >D); 

(3) , the samples are further cooled slowly (typically, < 10 
°C/min) along the liquidus to between -40 and -100 °C, a 
process called freeze concentration (D— >E); and (4), the 
samples are transferred into liquid nitrogen for long-term 
storage (E— >I). In this approach, the formation of extracellu- 
lar ice leads to freeze concentration of the unfrozen solutions 



by ejecting solutes and cells from the frozen to unfrozen 
phases. As a result, dehydration of cells in the unfrozen 
phase ensues, which minimizes intracellular water available 
for ice formation inside the cells so that the cells can enter 
the glassy phase (I) easily when transferring into liquid ni- 
trogen. This approach typically requires a specialized ma- 
chine usually called controlled rate freezer (CRF) to achieve 
freezing in a controllable manner. The time required for the 
slow-freezing process is typically hours. 

Vitrification by definition is ice free. In other words, no 
(or negligible) ice formation or freezing will occur in the 
sample during cooling [48,51,54]. Conventional vitrification 
(A— >F— >II in Fig. 3) has also been studied for both cells and 
tissue. In this approach, biological samples (state A) are first 
loaded with a very high concentration of CPA(s) (up to ~ 7 
M, state F) [48,51,54]. The samples are then cooled directly 
from ambient temperature to a cryogenic temperature usually 
in liquid nitrogen (state II) and stored there for future use. 
Although the conventional vitrification approach can be used 
to eliminate the detrimental effect of ice formation alto- 
gether, the unusually high concentration of CPA required by 
the approach is toxic to most mammalian cells even in a 
short period of exposure (ranging from seconds to minutes 
dependent on the specific cells and tissues) [55-60]. There- 
fore, the samples should be cooled as soon as possible after 
loading with CPAs. A mixture of multiple CPAs is often 
used to reduce the cytotoxicity of the high CPA concentra- 
tion required [61]. In addition, large, membrane imperme- 
able molecules such as sugars (typically sucrose and treha- 
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lose) have been used to minimize ice formation and protect 
cell membrane from injury during cooling [62-64]. Vitrifica- 
tion can be done without a specialized machine and the time 
required is generally much shorter than that for slow- 
freezing. 

Low-CPA vitrification (A— >B— ►III in Fig. 3) is a further 
advancement of the conventional vitrification with the goal 
to reduce the CPA concentration (e.g., at state F vs. B) re- 
quired for vitrification to a low, nontoxic level (similar to 
that used for slow-freezing). This can be done by creating an 
ultrafast cooling rate to cool the cells for cryopreservation. 
This is because the higher the cooling rate, the less the 
amount of cryoprotectants is required for achieving vitrifica- 
tion [64-69]. For example, even pure water can be vitrified 
without any ice formation when the cooling rate is about one 
million degree Celsius per second [69-71]. 

Various devices have been utilized to achieve fast cool- 
ing rates (~ 20,000 °C/min) such as the traditional French 
type straw, open pulled straw, electron grid, and cryoloops 
[56,72-74]. As a result of the fast cooling rate, the amount of 
cryoprotectants required for vitrification can be reduced to 
around 4 M. To achieve an ultrafast cooling rate, two recent 
studies reported the use of a micro-fabricated oscillating heat 
pipe (OHP) device [75,76]. Although their theoretical analy- 
sis shows that an ultrafast cooling rate of ~ 10 6 °C/min could 
be achieved, testing of the device for low-CPA vitrification 
using living cells has not been reported to date. Another re- 
cent study reported that a cooling rate as high as 200,000 
°C/min can be achieved by plunging an ultra- thin walled (10 
jam) quartz microcapillary (QMC, 180 jam inner diameter, 
slightly bigger than the diameter of a human oocyte) into 
liquid nitrogen [64]. As a result, the CPA concentration re- 
quired for vitrification of mouse embryonic stem cells and 
mouse oocytes was found to be as low as 2.5 M altogether 
[64,77], which is close to the upper boundary of CPA con- 
centrations used for slow- freezing. More recently, it was 
found that the required CPA concentration could be further 
reduced to ~ 1.5 M when mouse mesenchymal stem cells 
were encapsulated in ~ 100 urn alginate microcapsules for 
vitrifying using the thin-walled (380 jam inner diameter) 
QMC [78]. Therefore, the QMC-based low-CPA vitrification 
technique is promising to provide a solution to the major 
challenge facing the conventional vitrification approach as- 
sociated with cytotoxicity of high CPA concentration. Be- 
cause low-CPA vitrification avoids the shortcomings of both 
the conventional slow-freezing and vitrification approaches 
while combining their advantages, it is becoming increas- 
ingly popular, particularly for the cryopreservation of impor- 
tant and osmotically sensitive mammalian cells such as oo- 
cytes, sperm, stem cells, and many primary cells important 
for tissue engineering, cell-based therapy, and assisted re- 
production. 

1.2.2. Thermal Stabilization at Ambient Temperatures (Ly- 
opreservation) 

With the development of modern cell-based medicine, 
the demand on wide distribution of biopreserved products to 
end users in medical units at both urban and remote locations 
is increasingly high. The difficulty to widely distribute cryo- 
preserved biologicals in liquid nitrogen has motivated cryo- 
biologist and engineers to develop approaches for long-term 
stabilization of biologicals at ambient temperature. The idea 
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Fig. (4). Tardigrade (A and B) and resurrection plant (C and D) can 
survive extreme drought in nature upon rehydration - a phenome- 
non called anhydrobiosis or life without water. Figure reprinted and 
redrawn from references [85,86] with permission from John H. 
Crowe at U of California-Davis (for panels A and B) and Wayne P. 
Armstrong at Palomar College (for panels C and D). 

of dry or lyopreservation at ambient temperature is actually 
not new since many lower organisms (e.g., tardigrade shown 
in Fig. 4A and B), resurrection plants (e.g., Selaginella Lepi- 
dophylla shown in Fig. 4C and D), and seeds can survive 
extreme drought in nature upon rehydration, a phenomenon 
called anhydrobiosis or life without water [79-86]. A high 
concentration of sugars (typically sucrose for plants and tre- 
halose for lower organisms) has been found in these organ- 
isms and plants when they are in the anhydrobiotic state. 
Learning from nature, both sucrose and trehalose have been 
used as the protective agent (also termed lyoprotectant) in 
protocols for lyopreservation [87-91]. In practice, two ap- 
proaches have been studied to bring biologicals in aqueous 
samples to a dry state: freeze-drying (or lyophilization) and 
evaporative or convective drying (or desiccation) as shown 
in Fig. (3), as well. 

A typical freeze-drying protocol is as follows 
(A^B^C^D^G^H^IV in Fig. 3): (1), a biological 
sample such as proteins and cells in an aqueous solution 
(state A) is first supplemented with lyoprotectants (e.g. su- 
crose and trehalose) at a concentration of up to ~ 15 wt% 
(state B); (2), the sample is sub-cooled to usually between -3 
and -7 °C (phase C) to seed ice in the solution by touching 
the samples with a deeply cooled object (C— >D); (3), the 
sample is further cooled to between -30 and -50 °C slowly at 
a cooling rate usually less than 10 °C/min (D— >G); (4), the 
ice formed in the sample during freezing is then sublimated 
by exposing the sample to a vacuum usually less than 10 Pa 
at phase G (primary drying); and (5) a secondary drying 
process is then done by heating the sample in vacuum slowly 
to ambient temperature to further dehydrate the sample for 
additional hours to days (G— >H— >IV). The samples are then 
sealed and preserved in the dry phase (IV) at ambient tem- 
perature for future use. Freeze-drying has been used success- 
fully in achieving lyopreservation of many biomacro- 
molecules such as proteins and lipids and many pharmaceu- 
tical drugs. It has also been used for achieving lyopreserva- 
tion of prokaryotic cells such as bacterium, red blood cells, 
and platelets, but not eukaryotic cells at this time. For freeze- 
drying, it is crucial to keep the temperature low enough (be- 
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low the so-called collapse temperature [92-106]) during pri- 
mary drying (at phase G). Otherwise, the sample may col- 
lapse (can't maintain the morphology of the frozen sample) 
resulting in incomplete drying and heterogeneity in the 
freeze-dried product, which could significantly decrease the 
biostability of the freeze-dried biologicals [107-109]. More 
importantly, a recent study reported that the collapse tem- 
perature of cell culture medium-based trehalose solutions 
important for freeze-drying mammalian cells can be much 
lower than that of a simple binary trehalose-water solution 
and trehalose solutions used for freeze-drying pharmaceuti- 
cals and prokaryotes [110]. 

Unlike freeze-drying, during desiccation by evapora- 
tive/convective drying (A— >B— >IV in Fig. 3) water in an 
aqueous sample is removed by exposing the sample to a dry 
environment (e.g., dry air, inert gas such as nitrogen, and 
vacuum) without freezing (or ice formation) after loading 
with up to 15 wt% lyoprotectants. Desiccation by evapora- 
tive/convective drying has been used to achieve lyopreserva- 
tion of both biomacromolecules such as proteins and lipids, 
pharmaceutical drugs, and prokaryotic cells (such as bacte- 
rium, red blood cells, and platelets), but not eukaryotic 
mammalian cells. A major engineering difficulty to dry the 
glass-forming trehalose based solution for cell lyopreserva- 
tion by convective drying is that a thin glassy skin can easily 
form on the interface between the solution and the dry envi- 
ronment, resulting in incomplete drying and heterogeneity in 
the dried product. This problem might be minimized by 
breaking down the solution into micron- sized droplets or 
thin-films [108]. 

Beside the engineering challenge to effectively dry the 
trehalose solutions, effective delivery of the small hydro- 
philic lyoprotectants (e.g., trehalose and sucrose) into mam- 
malian cells has been challenging as the first step toward cell 
preservation at ambient temperature. This is because lyopro- 
tectant such as trehalose must be present both intra and ex- 
tracellular ly to provide the maximum protection during dry- 
ing, but mammalian cells lack a mechanism to synthesize 
trehalose endogenously and their plasma membrane is im- 
permeable to the sugars [111-113]. Over the past decades, a 
number of approaches have been explored to introduce treha- 
lose into living cells for preservation purpose. The most 
straightforward approach is to deliver exogenous trehalose 
into the cytosol of living cells by direct microinjection. This 
approach has been successfully used for intracellular deliv- 
ery of trehalose to cryopreserve mammalian oocytes that 
have a large size (-100 jam in diameter) and are generally 
in a small quantity (less than a few hundred) [113-116]. 

However, the microinjection approach is difficult (if not 
impossible) to apply for most living cells that are generally 
much smaller (< ~ 20 jam) than mammalian oocytes and usu- 
ally present in a large quantity (millions). Small living cells 
have been genetically engineered to synthesize trehalose 
endogenously. This approach requires the constant produc- 
tion of adenoviral vectors that exhibit significant cytotoxic- 
ity, particularly at high multiplicities of infection [117-119]. 
Trehalose has also been introduced into mammalian cells or 
their organelles through engineered or native transmembrane 
pores [91,112,120-122], electroporation [123,124], fluid- 
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phase endocytosis [125-127], and lipid phase transition 
[127,128]. 

In spite of the various approaches being explored, a con- 
sistent report of cell preservation using trehalose for small 
eukaryotic living cells is still absent [88,111,129,130]. This 
could be due to the inability to deliver a sufficient amount of 
intracellular trehalose (> ~ 0.1 M) for cell preservation using 
some of the approaches (e.g., fluid phase endocytosis). In 
addition, cells could be too severely compromised during the 
delivery step to withstand further freezing/dehydration 
stresses during preservation, considering the highly invasive 
nature of some of the approaches (e.g., electroporation). Re- 
cently, research has been sought to use liposomal and po- 
lymeric nanoparticles as the intracellular delivery vehicles of 
small hydrophilic molecules including the lyoprotectants 
with promising outcomes [131-135]. 

Besides the non-reducing disaccharide (trehalose and su- 
crose), small stress proteins particularly, the late embryo- 
genesis abundant (LEA) proteins have been suggested to be 
an important part of the molecular repertoire that renders 
desiccation tolerance in anhydrobiotic organisms and are 
attracting more and more research attention [79-84,136-141]. 

2. BIOLOGICAL EFFECT OF ABNORMAL TEM- 
PERATURES 

The biological responses to either abnormally low or 
high temperatures have been reported at the molecular, cellu- 
lar and tissue (both in vitro and in vivo) levels [142-149]. 
Some important structures/conformations of important bi- 




Fig. (5). An illustration (not to scale) of some important struc- 
tures/confirmations of a typical biological tissue (A), mammalian 
cell (B), cell plasma membrane (C), and protein (D): A biological 
tissue typically consists of cells (1), extracellular matrix (2) and 
blood vessels (3) through which oxygen and nutrients are trans- 
ported to the cells; a mammalian cell consists of many important 
sub-cellular organelles such as the rough endoplasmic reticulum 
(4), ribosome (5), smooth endoplasmic reticulum (6), mitochondria 
(7), plasma membrane (8), cytosol (9), lysosome (10), centriole 
(11), Golgi apparatus (12), and nucleus that are demonstrated as 
nuclear membrane (13), nucleolus (14), and nuclear membrane pore 
(15); the cell plasma membrane contains a phospholipid bilayer 
with hydrophilic heads (16) and hydrophobic tails (17), and mem- 
brane proteins (18) that works as transport channels, receptors, and 
structural connections between intracellular skeleton and extracel- 
lular matrix; and the most important secondary structures of a pro- 
tein are a helix (19) and (3 sheet (20). 
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Fig. (6). Fourier transform infrared (FTIR) and differential scanning calorimetry (DSC) signatures of thermally induced changes in cellular 
proteins in situ: (A), FTIR spectra at different temperatures showing the increase of extended (3 sheet structure and a simultaneous decrease 
of the a helix structure in cellular protein during heating mammalian cells slowly (2 °C/min); (B), DSC endotherm demonstrating heat ab- 
sorption as a result of protein denaturation during heating mammalian cells slowly also at 2 °C/min; (C), FTIR spectra showing extensive 
change of protein secondary structure from a helix to extended (3 sheet in thawed cells after freezing slowly to -80 °C, but not -20 °C; and 
(D), the DSC and FTIR measurements match well, suggesting protein denaturation is one of the major events that result in cell injury during 
heating between 40-70 °C. Figure reprinted and redrawn from references [158,195] with permission from Elsevier (for panels A and C) and 
Springer (for panels B and D). 



ologicals (e.g., tissue, cell, cell membrane, and protein) that 
have been reported to be the target of both cryo and hyper- 
thermic injury are shown in Fig. (5). The driving forces lead- 
ing to the biological responses at abnormally high tempera- 
ture are different from those at abnormally low temperatures 
due to the different biophysical events involved. The bio- 
logical responses at hyperthermic temperature appear to be 
dominantly thermotropic. At abnormally low temperatures, 
however, additional biophysical events such as ice formation 
and freeze concentration (see D— >E— >G in Fig. 3) as a result 
of ice formation at subzero temperatures can lead to addi- 
tional biological alterations at the molecular, cell, and tissue 
levels. 

2.1. Hyperthermic Injury 

A bevy of studies have been performed to investigate the 
alterations in biologicals exposing to hyperthermic tempera- 
tures. At the molecular level, lipid, protein, DNA, and RNA 
are the major structural and functional macromolecules of 
living cells and tissues. DNA and some structured RNAs do 
not undergo conformational change until above ~ 85-90 °C 
[150-153]. Therefore, they are unlikely to play an important 
role in cell killing at temperatures ranging from 37 to 85 °C 
for most thermal therapy protocols. Other RNAs that can 
perform enzymatic functions and form protein/RNA com- 
plexes such as ribosomes [151,154] could be potential tar- 
gets in thermal therapy applications [153]. Both protein de- 
naturation and lipid alteration (phase transition) have been 



commonly observed after thermal treatment [143]. Although 
many studies on protein denaturation were done using puri- 
fied proteins, the attention of recent studies have gradually 
shifted to study such changes in situ in living cells using 
various analytical techniques including differential scanning 
calorimetry (DSC), Fourier transform infrared (FTIR) spec- 
troscopy, and circular dichroism (CD) spectroscopy [155- 
159]. These studies have shown that alterations to lipid 
mainly occur below ~ 45 °C while major protein denatura- 
tion starts from 40 to 45 °C (Fig. 6 A and B) dependent on the 
heating rate and continue to occur at more than 100 °C. In 
view of the fact that major cell injury occurs at temperatures 
above 45 °C, protein denaturation has been proposed to play 
a more important role in defining cell injury than the change 
in lipids. This argument may be further strengthened by 
identifying the denaturation of specific proteins or protein 
groups that is the rate limiting step for thermal injury to cells 
[159]. However, lipid alterations are still considered to be 
important in defining the rate of hyperthermic injury because 
many cellular proteins (especially membrane proteins) are 
embedded in lipids and their stability is strongly affected by 
the lipids surrounding them [160-165]. 

Of note, with much attention being focused on studying 
the effect of hyperthermic temperatures on macromolecules, 
recent studies have found that a variety of molecules can be 
used to augment the performance of high temperatures for 
thermal destruction [166]. For example, the small peptide of 
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Fig. (7). Viability and function of mammalian cells post thermal treatment evaluated by exemplary assays: (A), immediate viability of cells 
assayed using calcein AM (green fluorescence indicating cellular metabolic activity and intact membrane of viable cells) and ethidium ho- 
modimer (red fluorescence indicating compromised cell membrane of injured and potentially dead cells); (B) and (C), expression of green 
fluorescent protein (GFP) under the control of Oct-4 gene (B) and a membrane surface glycoprotein (SSEA-1, C) indicating the undifferenti- 
ated properties of Rl murine embryonic stem cells post thermal treatment; and (D) histological difference between viable and damaged renal 
cells in normal porcine kidney tissue after 2 day's culture in media post thermal treatment showing the transition (arrow) from the lower left 
zone of significant necrotic tissue to the upper right zone with mainly intact tissue: Damaged cells lack nuclei. Figure reprinted and redrawn 
from references [38,64] with permission from Elsevier (for panels B and C) and Informa (for panel D). 



tissue necrosis factor- alpha (TNF-oc, 12.6 kD) has been 
shown to sensitize tumor for thermal therapy when delivered 
to the tumor site using functionalized gold nanoparticles 
[167]. Anticancer drugs encapsulated in liposomes have been 
shown to significantly increase the effect of hyperthermia in 
inhibiting tumor growth in vivo [168-170]. Using therapeutic 
adjuvants such as TNF-oc and anticancer drugs to augment 
thermal therapy and the targeted delivery of these adjuvants 
specifically into tumor using various nanoparticles (gold, 
magnetic, liposomal, and polymeric nanoparticles and car- 
bon nanotubes for heating, as well) are attracting more and 
more attention in the field of thermal therapy [43,170-176]. 

At the cellular level, thermal damage to essentially all the 
subcellular organelles (e.g., the plasma membrane, mito- 
chondria, Golgi apparatus, and endoplasmic reticulum in Fig. 
5B) have been reported in the hyperthermic literature, pre- 
sumably due to the damage to their two important building 
blocks, the cellular proteins and lipids as discussed above 
[142]. The plasma membrane that provides cells with the 
first protection from its extracellular environment has been 
the major focus of many studies on hyperthermic injury. Hy- 
perthermic injury to the plasma membrane usually is mani- 
fested as bleb formation and hyper-permeability to large 
molecules that are impermeable to the plasma membrane of 
intact cells (Fig. 7A). Damage to the plasma membrane and 
other sub-cellular organelles can further cause cells to lose 
their capability of attaching to a substrate and the subsequent 
proliferation and therefore, cell death ensues. Of note, when 
cells are exposed to a hyperthermic temperature below 43 
°C, they can adapt themselves to the sub-lethal temperature 



and become more thermally resistant to subsequent heat 
treatments, a phenomenon called thermotolerance [14,177]. 
Cells obtain this capability by increasing the concentration of 
a group of constitutive stress proteins called the heat shock 
proteins inside the cells [178-184]. Heat shock proteins are 
known to protect critical cellular proteins for normal cellular 
function under stress conditions including that as a result of 
abnormal hyperthermic temperatures [180,181]. 

At the in vitro tissue level, living cells are embedded in 
their native extracellular matrix (mainly collagen, Fig. 5A), 
which is different from the substrate or culture medium for 
attached and suspended cells, respectively. The effect of dif- 
ferent extracellular matrices on the thermal sensitivity of 
mammalian cells has been reported in a few studies [185- 
188]. The results however, are still inconsistent. Some stud- 
ies show that cells are more thermally sensitive in their na- 
tive matrix while others show no significant difference or 
even the opposite [185-188]. Further studies in a more con- 
trolled fashion are required to clarify the discrepancy in the 
existing literature in this regard. 

At the in vivo tissue level, blood perfusion (Fig. 5 A) and 
inflammatory wounding healing responses have been shown 
to play an important role in determining the extent of thermal 
injury: The thermal threshold required to destroy cells em- 
bedded in in vivo tissue is lower than that in vitro 
[38,167,189,190]. This observation has been attributed to 
vascular stasis as a result of thermal damage to tissue vascu- 
lature during and post thermal therapy, which leads to the 
deprivation of nutrition and oxygen, a secondary mechanism 
that kills the cells in addition to direct cell injury. Inflamma- 
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Fig. (8). Both histology and blood perfusion defects have been used 
to evaluate tissue injury in vivo: (A), a micrograph showing the 
patterns of in vivo hyperthermic injury in porcine renal tissue 7 
days after thermal treatment and (B) in vivo blood perfusion defects 
due to cryothermic injury in prostate tissue at different times (i.e. 0, 
3 and 7 days) measured in situ by fluorescence contrast. Although 
wound healing was observed in the necrotic and partially ablated 
zone, a central lesion named thermal fixation knowing for its resis- 
tance to wound healing response were observed in kidney tissue 
exposed to hyperthermic temperatures (A). Wound healing was 
observed throughout the cryogenic lesion leading to its shrinkage 
from day 3 to day 7 (B). Figure reprinted and redrawn from refer- 
ences [38,247] with permission from Informa (for panel A) and 
ASME (for panel B). 

tory wound healing response has been reported to be limited 
to the peripheral (i.e., the partial ablation and necrotic zone 
in Fig. 8 A) of the thermal lesion in vivo even after more than 
one week [38,191]. The lack of wound healing response in 
the central zone (i.e., the thermal fixation zone in Fig. 8A) of 
a hyperthermic lesion in vivo is presumably due to the lack 
of blood perfusion as a result of vascular stasis post thermal 
therapy in the observed time period. Sufficient blood perfu- 
sion is necessary for recruiting inflammatory cells such as 
neutrophils and macrophages for the wound healing response 
to occur [192]. As a result of the lack of wound healing re- 
sponse in the central lesion, it is reported that the host actu- 
ally takes the central lesion as a foreign body and tries to 
wall out the central lesion from its surrounding tissue by 
dystrophic calcification in the peripheral zone [38,191,193]. 
More studies are necessary to clarify the long-term (i.e., in 
months and years) fate of the central zone in a hyperthermic 
lesion in vivo. 

2.2. Cryothermic Injury 

The response of biologicals to cryogenic temperature 
may not be simply thermotropic in that two additional bio- 
physical events including ice formation and freeze concen- 
tration/dehydration may occur and can result in significant 
damage to cells [50,68,194]. As with hyperthermic injury, 
cryothermic injury has been investigated at the molecular, 
cellular and tissue (both in vitro and in vivo) levels [142- 
144]. At the molecular level, alterations in both lipid and 
proteins in response to cryogenic temperatures have been 
observed [195-197]. Dehydration and excursion in tempera- 
ture are the thermodynamic driving forces for lipid phase 
separation and phase change among various states including 
liquid crystal lamellar (hydrated), gel (less hydrated), and 
even inverted micellar (Hex II) in bulk lipid membranes 
(composed of many types of lipids) [144]. Both phase sepa- 
ration and the Hex II phase change in lipids have been re- 
ported to be highly correlated with cell membrane damage 
which ultimately destroys cells upon thawing in various 
plant systems [198,199]. Presumably, lipid phase change and 
separation are also important in defining cryothermic injury 
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to the plasma membranes of mammalian cells, although it 
has not been as clearly demonstrated as that in plant systems. 

Protein denaturation during freezing is also inevitable be- 
cause a freezing process is generally associated with change 
in pH and solute/electrolyte concentration, both of which can 
result in conformational alterations of proteins [143]. It has 
been reported that denaturation of cellular protein was mini- 
mal in thawed cells after freezing them to -20 °C (Fig. 6C) 
[195]. However, significant denaturation of cellular proteins 
was observed when the cells were frozen to -80 °C and 
thawed (Fig. 6C) [195]. In the latter case, an apparent shift 
of the secondary structure of cellular protein from a-helix to 
extended (3-sheet (Fig. 6C) is observable in the FTIR spectra, 
indicating denaturation and aggregation of cellular proteins. 
Interestingly, cell viability also correlates well with the de- 
naturation of cellular proteins at -80 °C according to the 
study [195]. Finally, as with hyperthermic injury, various 
molecules including TNF-oc [200-205], antifreeze proteins 
[206,207], excessive electrolytes [208,209], immunopoten- 
tiator [210-213], and various anticancer drugs [214-218] 
have been shown to augment the efficacy of low tempera- 
tures for thermal destruction. For example, TNF-oc has been 
found to increase the thermal threshold of cryoinjury to pros- 
tate cancer cells from ~ -20 to 0 °C [200], which is signifi- 
cant because it indicates that all cancer cells could be killed 
in a frozen tumor iceball and it is much more convenient to 
monitor the size of the frozen tumor iceball (e.g., using ultra- 
sonography) than the subzero temperature inside the iceball. 
Therefore, using therapeutic adjuvants such as TNF-oc to 
augment cryoinjury and the controlled/targeted delivery of 
the adjuvants into tumor and tumor cells using nanotechnol- 
ogy have attracted more and more attention in the field of 
thermal destruction at cryogenic temperatures [201,219,220]. 

At the cellular level, cryothermic injury can be attributed 
to ice formation induced mechanical damage and freeze con- 
centration induced physiochemical deviation from the 
physiological state. The latter can be connected to the altera- 
tion to cellular proteins and lipids at the molecular level as 
discussed previously. Ice formation can occur both intracel- 
lular ly and extracellularly. Although significant intracellular 
ice formation (IIF) is generally considered to be a lethal 
event, extracellular ice is not as damaging [50,194]. Signifi- 
cant IIF could destroy the cellular plasma membrane and 
other sub-cellular organelles mechanically as a result of 
volumetric expansion associated with the phase change or 
solidification of cellular water to ice, and the stochastic event 
of crack formation in the frozen cellular interior during 
freezing and thawing. The growing dendritic front of ice 
during IIF could also dismantle the native structure of cellu- 
lar proteins, lipids, and subcellular organelles. In addition, 
significant IIF could lead to freeze concentration of the cellu- 
lar interior which is detrimental to cellular proteins and lip- 
ids. Therefore, multiple mechanisms contribute to the lethal- 
ity of significant IIF to cells. 

For applications with a low concentration of CPA (e.g., 
cryopreservation by the conventional slow-freezing) or with- 
out CPA (e.g., thermal destruction or cryosurgery), signifi- 
cant IIF and freeze concentration induced excessive dehydra- 
tion are well-established to be the two primary biophysical 
events that cause cryoinjury. Since freeze concentration re- 
quires sufficient time for dehydration of cellular water across 
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the cell plasma membrane, its cell damaging effect is appar- 
ent only when cooling cells at slow cooling rates. For IIF, 
however, a higher cooling rate (tens to hundreds of degree 
Celsius per minute) is required so that intracellular water has 
no time to diffuse out of the cell and is available for IIF. This 
differential dependence on cooling rate of the damaging ef- 
fect of IIF and freeze concentration results in the classical 
inverted U curve (dark solid line (D in Fig. 9) of cell viability 
vs. cooling rate for the conventional slow-freezing (for cryo- 
preservation) and cryosurgery applications [50]. Of note, this 
inverted U curve does not apply to thermal stabilization by 
vitrification. For conventional vitrification, the cell injury 
mechanism is related to the high concentration of CPA used, 
which leads to cell damage by applying a significant osmotic 
stress and disturbing the normal metabolic pathways at su- 
per-zero temperatures before cooling and after warming the 
cells [55-60,221]. For low-CPA vitrification, cells loaded 
with a low concentration of CPA are brought to the vitrified 
phase without any ice formation or freeze concentration. 
Therefore, the traditional inverted U curve for cell survival 
vs. cooling rate can be modified and extended into the ul- 
trafast cooling rate domain where cell viability increases 
with increasing cooling rate (dark dashed line © in Fig. 9), 
which typically occurs when the cooling rate is as high as 
thousands to millions of degree Celsius per second when the 
CPA concentration is lower than 2-4 M. 

At the in vitro tissue level, the cell injury mechanism dis- 
cussed at the molecular and cellular level should still apply. 
In general, cells isolated from a tissue tend to transport water 
more quickly than cells within the intact tissue [144]. There- 
fore, cells form intracellular ice more readily in tissues due 
to enhanced water trapping as well as cell-ECM (extracellu- 
lar matrix) and cell-cell contact issues. For example, IIF in 
primary hepatocytes, fibroblast, smooth muscle cells were 



observed to be enhanced after placing them in a matrix (of 
either collagen or fibrin) over suspension [222-224]. IIF has 
been suggested to propagate through the gap junction be- 
tween adjacent cells [225-227]. Mechanically, abnormal 
thermal stress higher than the yield stress has been shown to 
accumulate in tissue and micro and macro crack formation in 
the tissue is commonly observed, which can lead to damage 
to the cells and the integrity of the extracellular matrix in 
both cryosurgery and cryopreservation applications [228- 
237]. Freezing has been shown to result in extensive damage 
to the vasculature in renal tissue, which presumably is re- 
sponsible for the side effect of significant bleeding post 
cryosurgical operation of the kidney [238-243]. Blood vessel 
is the target of freezing induced mechanical damage, pre- 
sumably because it resembles the defect in an otherwise 
more homogeneous tissue [229]. Freezing has also been 
shown to dramatically affect the biomechanical properties 
such as the porosity of the extracellular matrix of soft tissues 
[244] and the Young's modulus and ultimate tensile strength 
(UTS) of blood vessels [236]. Considering the tight connec- 
tion between cells and extracellular matrix through focal 
adhesion, mechanical damage to the extracellular matrix may 
further disrupt the integrity of cellular membrane. Of note, 
the presence of CPA can greatly reduce the mechanical dam- 
age, particularly when the CPA concentration is high enough 
to result in vitrification [245,246]. 

At the in vivo tissue level, several studies have shown 
that the subzero temperature required for causing tissue ne- 
crosis is higher than that in vitro [200,239,247]. Again, this 
observation has been attributed to the vascular stasis as a 
result of freezing induced damage to tissue vasculature dur- 
ing freezing and thawing, which results in the deprivation of 
nutrition and oxygen to cells, a secondary mechanism that 
kills cells in addition to direct cell injury. The wound healing 



Thermostability of Biological Systems 

response in an in vivo cryolesion however, is quite different 
from that in an in vivo hyperthermic lesion in that wound 
healing has been observed throughout the whole cryolesion 
[191]. For example, the cryolesion shrinks from day 3 to day 
7 (Fig. 8B) indicating wound healing response that assimi- 
lates the necrotic tissue and allow the regeneration of normal 
tissue in the cryolesion [144,247]. This is probably because 
the frozen tissue has generally been observed to be reper- 
fused immediately after thawing and vascular stasis does not 
occur usually until a few hours after thawing in cryosurgical 
applications [239,247]. Presumably, the short period of tem- 
porary reperfusion of the frozen tissue after thawing allows 
the recruitment of inflammatory cells for the subsequent 
wound healing response in the whole damaged tissue do- 
main. 

2.3. Measurement of Thermal (Cryo and Hyperthermic) 
Injury 

Various methods have been utilized to measure thermal 
injury to biologicals. At the molecular level, x-ray crystal- 
lography, nuclear magnetic resonance (NMR) spectroscopy, 
Fourier transform infrared (FTIR) spectroscopy, circular 
dichroism (CD), electron microscopy, and differential scan- 
ning calorimetry (DSC), have been used to investigate ther- 
mally induced alterations to many proteins and lipids and 
sometimes DNA/RNAs [143,187]. Both x-ray crystallogra- 
phy and NMR spectroscopy can determine definitively the 
3D structure of pure proteins. However, x-ray crystallogra- 
phy requires the protein sample to be in the crystal state and 
NMR spectroscopy can only be used to study small proteins 
or protein domains (< 20 kD) [151]. Both FTIR (Fig. 6A and 
C) and CD have been widely used to study the structural 
change of proteins (secondary structures such as a-helix and 
(3-sheet in Fig. 5D) and lipids in response to a thermal chal- 
lenge [143,187]. DSC is the only technique that can be used 
to investigate the calorimetric effect associated with the 
structural change in proteins (denaturation and agglomera- 
tion, Fig. 6B), lipids (phase transition and separation), and 
DNA/RNAs (denaturation, breaking up the intra- molecular 
hydrogen bond). Although isolated proteins and lipids have 
been used in many studies of thermal injury, there is a shift 
of research interests to investigate the in situ protein denatu- 
ration and lipid change in living cells recently. This is impor- 
tant because protein and lipid stability is strongly affected by 
the biochemical microenvironment in which they reside 
[143]. This is particularly important for determining the cor- 
relation between cell injury and thermally induced alterations 
in proteins and lipids. A number of studies using FTIR and 
DSC in this regard have reported strong correlations between 
cell injury and in situ thermal denaturation of protein in a 
number of cell types including prostate cancer cells (Fig. 
6D), primary liver cells, and red blood cells 
[156,158,159,195]. 

Immediate cell viability (or injury) post thermal treat- 
ment is usually judged by the integrity of cell plasma mem- 
brane. The ability of an intact plasma membrane to exclude 
some membrane impermeable dye (e.g., ethidium homodi- 
mer showing red fluorescence in Fig. (7A), propidium io- 
dide, and trypan blue) has been the most convenient way to 
visualize thermally induced damage to mammalian cells us- 
ing microscopy (fluorescence capability is required for using 
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ethidium homodimer and propidium iodide) [187,248]. The 
hydrolysis of calcein-AM (a non- fluorescent hydrophobic 
compound that easily permeates intact cells) by intracellular 
esterases (important enzymes for intracellular hydrolytic 
reactions) to produce calcein (a hydrophilic compound with 
strong green fluorescence that is impermeable to the plasma 
membrane of intact cells) has also been widely used to judge 
immediate cell viability using fluorescence microscopy (Fig. 
7A). However, it has been shown that many cells with com- 
promised membrane (judged by impermeable dyes) also 
show strong green fluorescence of calcein. Therefore, cau- 
tion should be taken when calcein AM is used solely to 
judge immediate cell viability [187]. 

Since thermal damage to the other subcellular organelles 
can also determine the eventual fate of cells exposed to ab- 
normally high and low temperatures, a better assay to evalu- 
ate the true (or long-term) cell viability can be done by 
checking the capability of the cells to attach to a substrate 
and subsequently proliferate [187]. The former is usually 
termed attachment efficiency and is evaluated by calculating 
the percentage of attached cells with regard to total cells 
seeded 3-12 hrs after thermal treatment and seeding. Prolif- 
eration is obtained by monitoring the change in cell number 
in a sample in several consecutive days after treatment. 
Clonogenics, another assay for determining the long-term 
cell survival, is calculated as the ratio (or percentage) of the 
number of cell colonies formed for thermally treated cells to 
that of control (i.e., no thermal treatment) samples usually 
after 7-10 days of culture after seeding [249]. This assay 
avoids counting cells in multiple days for the proliferation 
assay. For thermal stabilization applications, it is also impor- 
tant to verify that the cells retain important cell specific func- 
tions. For example, it is important to make sure hepatocytes 
(or liver cells) retain the normal function of urea production 
and albumin synthesis after biopreservation [222]. For stem 
cells, it is important to make sure that they retain pluripo- 
tency (e.g., the undifferentiated properties indicated by the 
expression of some stem cells specific markers shown in Fig. 
7B and C for murine embryonic stem cells), multipotency, 
and/or the capability to differentiate into more specialized 
adult cells post biopreservation [250]. 

Histology (hemotoxylin&eosin or H&E stain) is com- 
monly used to judge cell (Fig. 7D) and vascular injury in 
tissue [192]. Although being qualitative initially, histology 
has been used to quantify thermal injury to cells in tissue in 
several recent studies (Fig. 8A) [38,186,251]. Membrane 
impermeable fluorescence dyes have also been used to quan- 
tify cell injury in tissue after thermal treatment [186,251]. 
For some tissues abundant of collagen such as rat tail tendon, 
joint capsule, and myocardium, thermal injury is reflected by 
collagen denaturation which can be measured by a variety of 
techniques such as DSC, loss of birefringence observable 
using polarized light microscopy, loss of transparency, and 
shrinkage as reviewed elsewhere [143,148,252]. For injury 
assessment in situ in the host, perfusion defect (or absence of 
blood perfusion) in the damaged tissue is commonly used, 
which can be measured by MR contrast agent (e.g., gadolin- 
ium) exclusion or by fluorescence contrast in situ (Fig. 8B) 
[200,253-255]. Recent studies have shown that tissue injury 
assessed by perfusion defect corresponds quantitatively with 
tissue necrosis determined by histology using biopsy 
[256,257]. 
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Fig. (10). An illustration of energy-state relationship of the 1 st order 
kinetic process of thermally induced protein denaturation: The na- 
tive protein is activated after absorbing energy followed by denatu- 
ration and potentially aggregation and agglomeration: Aggregation 
and agglomeration make the process irreversible. Figure reprinted 
and redrawn from reference [143] with permission from John Wiley 
& Sons. 

Among the various methods, clonogenics (or prolifera- 
tion) and histology are the gold standard for judging cell and 
tissue viability, respectively. Besides the methods used to 
determine it, the cell viability (0-100%) is also dependent on 
many other factors including the specific thermal history 
applied, the attachment state (i.e., suspended, attached on 
dish, attached in 3D matrices) of the cells, and the cell cul- 
ture microenvironment (e.g., pH) [187]. 

3. BIOPHYSICAL MODELING OF THERMAL IN- 
JURY 

Besides injury measurement, there is a need of injury 
models to predict the outcome of thermal treatment for better 
understanding of the injury mechanism in research and pre- 
treatment planning in the clinic. Injury models need to be 
verified or parameters in the model need to be extracted us- 
ing experimental data a priori. With the verified model, in- 
jury to biologicals under any thermal history can be pre- 
dicted. Furthermore, models may help us understand the in- 
jury mechanism by comparing the model parameters with 
those of known biological processes such as protein denatu- 
ration, lipid alteration, and DNA breakage [142,146,252]. 
Modeling hyperthermic injury at the molecular, cellular and 
tissue level has been reported in the literature for a wealth of 
biologicals [142,146,252]. However, modeling of cry other - 
mic injury is not as well-developed due to the complication 
of two additional biophysical events: freeze concentration 
(leading to cell dehydration) and IIF in addition to the ther- 
motropic effect of cryogenic temperature per se. 

3.1. Modeling of Hyperthermic Injury 

Several models have been developed to predict hyper- 
thermic injury to biologicals including the statistical model, 
enzyme denaturation model, thermal isoeffective dose (TID) 
model, and various kinetic models [142,146,258-261]. 
Among the various models, the TID and kinetic models are 
the most widely used. 



3.1.1. Kinetic Models 

The first order irreversible kinetic model has been the 
most commonly used for predicting hyperthermic injury. In 
this model, biologicals including biomacromolecules, cells, 
and tissue are considered to be either in a native (N) or in- 
jured (I) state with an injury rate (or rate constant in the 
chemistry literature) k as follows: 



N- 



^1 



(1) 

An energy- state illustration of the first order kinetic 
model for describing the denaturation of proteins is given in 
Fig. (10), where the denatured state of proteins is equivalent 
to the injured state of cells or tissue in Equation 1 . Thermally 
induced protein denaturation is often irreversible as a result 
of aggregation/coagulation and agglomeration following 
denaturation, which eventually results in irreversible cell 
death. 

Assuming all biologicals being in the native state ini- 
tially, the fraction of biologicals in the native state after be- 
ing exposed to a thermal history T(t) can be calculated as 
follows : 



F N = exp(-£fa/f) 



(2) 



where F represents fraction, t is time, and the subscript 'N' 
represents native state. The injury rate (or rate constant) can 
be calculated using the Arrhenius model as follows [262- 
264]: 



k = Atxp[-AE/(R g T)] 



(3) 



where R g is the universal gas constant (8.314 J mole" 1 K" 1 ), T 
(in K) is the thermal history (time dependent), and AE and A 
are the two kinetic model parameters (constants) usually 
called activation energy and frequency factor, respectively. 
This equation can be rewritten in the following format by 
taking the natural logarithm on both sides of the equation: 



\n(k) = \n(A)-AE/(R g T) 



(4) 



Equation 4 shows that the injury rate (k) is constant for a 
given temperature. Therefore, Equation 2 can be rewritten as 
follows under isothermal heating conditions: 



F N = exp(-fa) 



(5) 



Equation 5 has been frequently used to fit to data of frac- 
tional cell survival vs. isothermal heating time under various 
constant temperatures, through which the cell injury rate at 
various constant temperatures can be determined. With the 
known injury rate at various constant temperatures, the two 
model parameters (i.e., AE and A) can be determined by fit- 
ting to the data of k vs. T 1 at various constant temperatures 
using the linear relationship between the natural logarithm of 
k and T 1 given by Equation 4. The frequency factor and ac- 
tivation energy can be determined from the intercept and the 
slope of the linear fit, respectively. The procedures for de- 
termining the activation energy and frequency factor are 
demonstrated in Fig. (11). 

In the hyperthermic literature, a terminology called injury 
and often represented by Q is defined as follows [262-264] : 
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Fig. (11). An illustration of the procedures for determining the activation energy and frequency factor using Equations 4 and 5: (A), deter- 
mining the injury rate at various temperatures by fitting Equation 5 to experimental data on cell survival vs. exposure time and (B) determin- 
ing activation energy and frequency factor by fitting Equation 4 to the natural logarithm of injury rate vs. the inverse of temperature obtained 
from (A). Figure reprinted and redrawn from reference [187] with permission from ASME. 



Q= fkdt 
Jo 



(6) 



According to Equation 6, k is the first derivative of injury 
(Q) with respect to time, which explains why k is often 
called injury rate in the hyperthermic literature. The relation- 
ship between Q and the fraction of biologicals in the native 
state for a thermal treatment can be determined by combin- 
ing Equations 6 and 2 as follows: 



F N = exp(-Q) 



(7) 



From this equation, it is clear that an accumulated injury 
of 1 (i.e., Q=l) implies that the fractional survival is e" 1 or ~ 
36.7%. 

The reaction rate k can also be calculated using the abso- 
lute rate theory model as follows [265,266]: 



k = (k'/h')Tex V [-Ag"/(R g T)] 
Ag" = Ah" -TAs" 



(8) 



(9) 



where k' is the Boltzmann constant (1.38xl0" 23 J K" 1 ), h' is the 
Planck's constant (6.625xl0 34 J s), and Ah", As", and Ag" 
are enthalpy, entropy, and Gibb's free energy of activation 
for the thermal injury process, respectively. The two kinetic 
model parameters for the absolute rate theory model are the 
enthalpy and entropy of activation. The relationship between 
the parameters of the absolute rate theory model and those of 
the Arrhenius model for a first order reaction such as that 
described by Equation 1 is as follows [148]: 



AE = Ah" + R g T 
A = (k'/h')Texp(As"/R g +1) 



(10) 



(11) 



According to Equation 10, the difference (i.e. R g T) be- 
tween AE and Ah" is less than 3 kJ/mole when temperature is 
below 80°C. This small difference is negligible for most 
situations due to the relatively high value of the activation 
energy for cell/tissue injury and protein denaturation (typi- 
cally, more than 50 kJ/mole) [142]. The frequency factor is 
linearly dependent on temperature and exponentially de- 
pendent on the entropy of activation. Since the unit for tem- 



perature is Kelvin, it is usually more than 293.15 K (assum- 
ing room temperature is 20 °C) for thermal therapy applica- 
tions. Therefore, the dependence of the frequency factor on 
temperature is weak considering that the temperature for cell 
killing is usually between 50 and 70 °C in most thermal ther- 
apy applications. 

Many studies have been performed to determine the pa- 
rameters in the first order kinetic model for many different 
types of biologicals at the molecular (mainly proteins), cell, 
and tissue levels. These kinetic parameters for a wealth of 
biologicals have been summarized/tabulated and the follow- 
ing conclusion can be drawn by a careful examination of the 
kinetic data (Fig. 12) [142]: 




400 



800 1200 
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Fig. (12). A summary of the kinetic parameters in the first order 
kinetic model for various proteins, cells, and tissues: A simple, 
phenomenological, linear relationship exists between activation 
energy and the natural logarithm of frequency factor in the first 
order kinetic model for hyperthermic injury. Figure reprinted and 
redrawn from reference [142] with permission from Begell House. 

• The literature values of activation energy varies over a 
broad range from 50 to 1600 kJ/mole dependent on the type 
of biologicals (i.e., lipids, proteins, cells and/or tissues), 
temperature range, pH, and assay used. The activation en- 
ergy for lipid change is typically lower than that for protein 
denaturation and cell/tissue injury. The activation energy for 
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protein denaturation and cell/tissue injury however, is on the 
same order of magnitude, suggesting protein denaturation 
should be an important mechanism of thermally induced 
cell/tissue injury. 

• Clonogenic assay for single cells usually results in very 
high activation energy (~ 350-1600 kJ/mole) whereas it is 
generally much lower for histology and fluorescence viabil- 
ity dye assays (~ 100-800 kJ/mole). 

• An approximate, phenomenological, linear relationship 
between the natural logarithm of frequency factor and activa- 
tion energy exists (see Fig. 12). This approximate linear rela- 
tionship is attributed to the small variation of Gibbs free en- 
ergy of activation given by Equation 9 for a thermal injury 
process despite the large variation of activation enthalpy 
(related to activation energy by Equation 10) and entropy 
(related to the frequency factor by Equation 11). The Gibb's 
free energy of activation (Ag") typically varies between 100 
and 115 kJ/mole. 

• Temperature break points below and above which the 
activation energy changes significantly are mainly located at 
~ 43°C and ~ 50°C if they do exist [142,187]. The break 
point at 43 °C has been suggested to be a result of thermo tol- 
erance. The reason for a break point at 50°C is unclear yet. 

• With few exceptions, the literature values of frequency 
factor are greater than 2.0 x 10 13 s _1 which is larger than the 
frequency factor predicted using Equation 11 assuming 0 
activation entropy. In other words, the activation entropy for 
the kinetic process of thermally induced alterations in pro- 
tein, cells, and tissue are typically positive. 

Besides the first order kinetic model, higher order kinetic 
models have also been developed for modeling thermal in- 
jury [259,260]. However, higher order kinetic models have 
not been widely used probably due to the difficulty to defini- 
tively determine the multiple (> 2) parameters in those mod- 
els using experimental data. 

3.1.2. Thermal Isoeffective Dose (TID) Model 

The thermal isoeffective dose (TID) model was proposed 
nearly 30 years ago [267] and has been widely used in both 
the traditional and many recent hyperthermic literatures 
[149]. In the TID model, a thermal history (i.e., time- 
temperature data) is converted to a cumulative equivalent 
minutes (isoeffect) at 43 °C. As a result, the comparison of 
thermal doses and treatment planning of thermal therapy 
become straightforward, which may explain why the TID 
model has been widely used [268,269]. The time required to 
achieve the isoeffect of 43 °C at a constant temperature in 
terms of thermal injury is calculated as follows [149]: 

t(6) = CEM43C s /R 43 -° (12) 

where 6 is temperature in °C, CEM43C S is the cumulative 
equivalent minutes required to decrease cell survival from 
the initial (100%) to a given criterion (S can be calculated as 
F N x 100%) at 43 °C, and R is the number of minutes re- 
quired to compensate for a 1 °C temperature change to 
achieve an isoeffect of cell injury. The parameter R (not di- 
mensionless) is a function of temperature and the activation 
energy of a thermal injury process [267] : 

R = exp{-AE/[R g T(T + l)]} (13) 



If assuming 36.7% (e 1 ) cell survival as the criterion, the 
threshold time CEM43C can be calculated from the Ar- 
rhenius model as follows: 

CEM43C 361% =l/k 43C =A~ l exp{AE/[/^(43 + 273.15)]} 

(14) 

Note that CEM43C 36 .7% is the same as the inverse of the 
injury rate at 43 °C (k 43C ). To calculate CEM43C assuming a 
different cell survival S (CEM43C S ), the following equation 
can be used [187]: 

CEM 43C S = CEM43C 361% ln(S~ l ) (15) 

where the subscript S and 36.7% represent the percentage of 
survival that is assumed as the criterion for calculating the 
cumulative equivalent time at 43 °C. 

In many traditional literature of hyperthermia, a simpli- 
fied version of the TID model was used that assumes a uni- 
versal CEM43C of either 120 or 240 min for complete injury 
of all cells at 43 °C. The parameter R in the simplified TID 
model is taken as 0.25 and 0.5 for temperatures below and 
above 43 °C, respectively. According to Equations 13-15, 
both the CEM43C and R are dependent on the activation 
energy and frequency factor (R is also temperature depend- 
ent), both of which are strongly cell type dependent. There- 
fore, caution should be taken when using the simplified TID 
model for predicting thermal injury in different types of 
cells. The CEM43C for a variety of cells and tissue of 
urologic origin together with the temperature and cell type- 
dependent R values for thermal therapy applications have 
been tabulated elsewhere [187]. 

3.2. Modeling of Cryothermic Injury 

Modeling of cryothermic injury is more challenging since 
it is complicated by two additional but important biophysical 
events: freeze concentration and ice formation. To the 
author's knowledge, quantitative and mechanistic modeling 
of the effect of low temperature, freeze concentration, and 
ice formation on the stability of biomacromolecules such as 
lipids and proteins has not been reported. At the cell and 
tissue levels, formation of extracellular ice leads to freeze 
concentration along the liquidus in the phase diagram (see 
Fig. 3) which further results in cell dehydration at a slow 
cooling rate (e.g., < 10 °C/min). On the other hand, IIF will 
occur when the cooling rate is fast so that the intracellular 
water has no time to leave the cells through exosmosis, but 
not fast enough (< CR V shown in Fig. 9) to supersede the 
kinetics of ice formation. As with the molecular level injury, 
quantitative and mechanistic modeling of cryoinjury at the 
cellular and tissue levels has not been well studied. So far, 
much of the effort in modeling cryoinjury has been focused 
on modeling intracellular ice formation (IIF) and cell dehy- 
dration during freezing. IIF has been studied using both phe- 
nomenological and mechanistic models [194,270]. The 
mechanistic model has been widely used and delineates the 
ice formation as two consecutive events: Nucleation of ice 
nuclei and the subsequent growth of the nuclei [194,271]. 
Nucleation of intracellular ice can be catalyzed by either a 
surface (surface catalyzed nucleation, SCN) such as the cell 
plasma membrane or a volume of subcooled solution (vol- 
ume catalyzed nucleation, VCN) such as the cytoplasm 
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[194,271]. The rate of ice nucleation (I) due to either VCN or 
SCN can be modeled as follows [68,272]: 



I = Q f 
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where T f is the equilibrium freezing point of the intracellular 
solution, N is the number of water molecules either in the 
cells (for VCN) or in contact with the cell plasma membrane 
(for SCN), r| is viscosity, and Q and k are two model pa- 
rameters that are usually called the kinetic and thermody- 
namic coefficient, respectively, and the subscript 0 repre- 
sents the isotonic solution state before freezing (ice forma- 
tion). The two model parameters under isotonic solution 
state (Q 0 and K 0 ) need to be determined by experimental 
studies and have been reported in the literature for a number 
of cell types as summarized elsewhere [142,194,273]. The 
cumulative probability of intracellular ice formation (Puf) 
can then be calculated as follows [68,194]: 



P™=l-exp|-| VI™ dt 
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(18) 



(19) 



where V and A are the cell volume and surface area available 
for catalyzing the nucleation of intracellular ice, respectively. 
Significant IIF is usually manifested as darkening of the cell 
cytoplasm when observed under a bright field of light mi- 
croscopy as demonstrated in Fig. (13), which has been used 
widely to quantify the kinetics of intracellular ice formation 
[274,275]. 

Of note, the above IIF model is valid only when ice nu- 
cleation (the first step of IIF) is the rate- limiting step of IIF 
which is often true when freezing cells/tissues at not very 
high cooling rates (e.g., less than a few hundred Celsius per 
minute) and in the absence of a high concentration of CPA 
(e.g., less than 10 wt%) [194]. Under these conditions, ice 
will immediately propagate (or grow) throughout the entire 
intracellular space once the ice nuclei form, which results in 
the darkening event. When cooling cells/tissue with much 
higher cooling rates and/or a high concentration of intracel- 
lular CPA (e.g., during vitrification and after significant cell 
dehydration during slow-freezing), the rate limiting- step of 
IIF is the growth of the ice nuclei (the second step of IIF). 
The IIF under these conditions is said to be diffusion-limited 
for which more complicated models are needed to predict the 
amount of ice in cells [67,69,276-278]. To accurately predict 
the diffusion-limited ice nucleation and growth, an advanced 
model such as the free volume model that can be used to 
predict the effect of glass transition on solution viscosity and 
diffusion coefficient might be necessary. Such free volume 
model for a couple of cryo and lyoprotectants has been re- 
ported in [107]. 

In order to predict the probability of IIF using the above 
model during slow-freezing where freeze concentration in- 
duced cell dehydration is significant, information on the cell 
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Fig. (13). Dependence on cooling rate of intracellular ice formation 
(IIF) when freezing hepatocytes from -4 (after seeding extracellular 
ice) to -40 °C showing the darkening of the cellular interior as a 
result of significant IIF under a bright field of light microscope: 
When the cooling rate is very low (10 °C/min), cell dehydration 
dominates (A); when the cooling rate is high (400 °C/min), IIF 
dominates (C); and when the cooling rate is intermediate (100 
°C/min), both IIF and cell dehydration can occur (B). Figure re- 
printed and redrawn from reference [272] with permission from 
John Wiley & Sons. 

volume, V (or cell surface area, A, related to the diameter of 
the cells when the cells assume a spherical geometry), during 
freezing is required. The following model has been com- 
monly used to predict the cell volume change during freez- 
ing [67,279]: 



dV_ 
dt 



L p AR g T 
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(20) 



where v w is the partial molar volume of water, n represents 
amount (in mole), (p s is dissociation constant of solutes (e.g., 
2 for NaCl), V s is volume of solutes, V b is osmotically inac- 
tive volume in cells, Ah f is latent heat of fusion of water, T ref 
is a reference temperature usually taken as the equilibrium 
melting point of the intracellular solution, T is thermal his- 
tory, the subscripts s and w represent solute (including 
CPAs) and water, respectively, and L p is the cell plasma 
membrane permeability to water that can be calculated as 
follows [280]: 
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(21) 



where L pg is the permeability of the cell membrane to water 
at the reference temperature (T re f) and E Lp is the activation 
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energy for the water transport across the cell plasma mem- 
brane. In the equation, L pg and E Lp are two model parameters 
that need to be determined a priori using experimental data. 
Cell dehydration during freezing can be measured using ei- 
ther a specialized cryostage mounted on a light microscope 
[275] or differential scanning calorimetry [281-284]. Many 
studies have been performed to determine the two model 
parameters for various cells, which were reviewed and tabu- 
lated elsewhere [142,273,284,285]. 

Although cell dehydration and intracellular ice formation 
(IIF) can be predicted using the above models, a quantitative 
understanding of the mechanistic link between the two bio- 
physical events and cell injury has not been well established 
despite some early efforts in this respect as reviewed else- 
where [142]. The incidence of significant IIF (defined as 
darkening of cell cytoplasm) correlates strongly with cell 
death in many cell types (i.e., 50% of IIF in many cell popu- 
lations yields 50% of dead cells) [194]. However, the exact 
amount or percentage of intracellular ice that is significant 
enough to result in irreversible cell death is still unclear. 
Some studies even suggest that a small amount of intracellu- 
lar ice might be beneficial to cell survival [286,287]. There- 
fore, further studies to establish mechanistic models capable 
of accounting for the effect of all the freezing induced bio- 
physical events including intracellular ice formation and 
freeze concentration (i.e., the so-called solute effect), and 
low temperatures per se is important to further our under- 
standing of low temperature biology and its biomedical ap- 
plication such as cryosurgery and cryopreservation. 

4. THERMOPHYSICAL QUANTIFICATION 

Heat transfer in living systems is unique because of the 
effect on heat transfer of blood perfusion (and sometimes 
metabolism) (Fig. 5 A) [21]. Generally, the thermal response 
in a bioheat transfer process (either heating or cooling) can 
be modeled using the following energy balance equation 
together with appropriate initial and boundary conditions: 



d(cT + AL) 
' dt 



= V.(k T VT) + Q EE +Q BP +Q M (22) 



where p is density, c is specific heat, L is latent heat (per unit 
mass), a is the latent heat release pattern that could be tem- 
perature and/or time dependent [142], k T is thermal conduc- 
tivity, Q represents heat source/sink term, and the subscripts 
EE, BP, and M represent the source/sink term due to the ap- 
plied external energy, blood perfusion, and metabolism, re- 
spectively. A detailed discussion about the different terms is 
given below. Note that the effect on heat transfer due to con- 
vection is omitted in the above energy equation since con- 
vection in the bulk (excluding the blood vessel that is ac- 
counted for by the heat source/sink term due to blood perfu- 
sion (Qbp)) soft tissue is generally considered to be negligi- 
ble. 

4.1. Metabolic Heat Generation 

In many studies of thermal therapy and cryosurgery, 
metabolic heat generation is taken as zero because it is gen- 
erally much smaller than that due to blood perfusion and the 
externally applied energy source. In some studies particu- 
larly for the modeling of heating or cooling without freezing 
(no ice formation), the metabolic source term is taken to be 
temperature dependent as follows [288,289]: 



Q M =Q M ^ (mo (23) 

where Qm,37c is the baseline metabolic heat generation at 37 
°C that is tissue type dependent and is summarized in [146] 
for a number of tissues/organs, $ is a model parameter that 
should be tissue type dependent as well although it is usually 
taken as 3 in the literature [288,289]. 

4.2. External Energy Source Term 

The external energy source term (Qee) is related to the so- 
called specific absorption rate (SAR) in the literature, par- 
ticularly when RF/microwave is used as the energy source 
for thermally destroying unwanted tissue or warming cryo- 
preserved biologicals. The unit of Q E e is power per volume 
while SAR is in power per weight. Therefore, Qee is the 
product of SAR and density. The direct approach for deter- 
mining this external source term is to solve either the elec- 
tromagnetic (for microwave, RF, and laser) or mechanical 
(for high intensity focused ultrasound or HIFU) wave equa- 
tion to determine the electrical or pressure field (i.e., E and 
p), respectively [35,290,291]. For RF/microwave applica- 
tions, Qee can be calculated from the electric field as follows 
[40]: 



Q ee =(o + (De")\e\ 



(24) 



where a is electrical conductivity, co is angular frequency, e" 
is the imaginary part of electrical permittivity, and || E || rep- 
resents the magnitude of electric field. For laser and HIFU 
heating, the governing equations for wave propagation be- 
come very nonlinear as a result of wave scattering and cavi- 
tation (for HIFU). Therefore, a simplified model derived 
based on the well-known Beer-Lambert law of absorption 
has been utilized to predict the Q E e of a plane wave as fol- 
lows [292]: 



Q EE =aI 0 e- az f(r) 



(25) 



where a is absorption coefficient that accounts for all the 
mechanisms of heat generation, I 0 is the incident wave inten- 
sity, z is the coordinate in the direction of laser or HIFU 
wave propagation, and f(r) is a function of the intensity dis- 
tribution in the radial direction which is usually taken in the 
form of Gaussian distribution as follows [292] : 



/(r) = exp[-2(r/r 0 ) 2 ] 



(26) 



where r is the radial coordinate and r 0 is known as the 'e~ 
radius' of the wave beam since at r = r 0 , f(r)=e" 2 . 

Besides theoretical prediction, experimental measure- 
ment of SAR (or Qee) has also been widely adopted for 
complicated probe designs, particularly when RF/microwave 
is used as the energy source. Since the experimental method 
can only determine SAR values at discrete, representative 
locations, analytical expressions are usually used to fit to the 
experimental data to predict SAR in a continuous domain 
[38,39]. For example, Fig. (14) shows the measured SAR 
pattern together with its fit for a typical axisymmetric dipolar 
helical microwave probes using the following equations 
[38]: 
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Fig. (14). A typical SAR pattern of a dipole microwave probe showing the penetration of microwave energy into the tissue space (in r direc- 
tion). The frequency of the microwave used is 921 MHz. Figure reprinted and redrawn from reference [38] with permission from Informa. 



SAR = Q E /p = Ae- 2(5r r- n (l + a \z\ + bz 2 
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(27) 



where r is radial coordinate originated from the probe center- 
line, z is the axial coordinate originated from the middle of 
the microwave thermal probe, (3 (= 0.413 cm" 1 for kidney 
tissue) is attenuation constant, |z| represents the absolute of z, 
and the constant model parameters A, n, a, b, c, d, and z 0 are 
1.24xl0 3 , 0.9523, -0.3065, 2.783, -0.941, 0.2295, and 
1.0515, respectively. The two SAR peaks at r = 0 are the 
result of the dipole design of the microwave thermal probe. 
The SAR distribution along the axial direction of the probe 
(except near the two ends) is more homogeneous than that 
along the radial direction. SAR generally decreases exponen- 
tially in the radial direction due to tissue absorption of the 
electromagnetic energy. 

4.3. Thermal Effect of Blood Perfusion 

The effect of blood perfusion on thermal modeling in un- 
frozen tissue particularly during heating has been the focus 
of many studies in the literature. In frozen tissue however, 
blood perfusion is stopped and has no thermal effect. The 
classical and probably the most frequently used model of the 
blood perfusion term is given as follows (note, Equation 22 
with the blood perfusion term given below is often called the 
classical Pennes bioheat equation) [293]: 



Q BP =wp B c B (T a -T) 



(28) 



where w is blood perfusion rate and the subscripts B and a 
represent blood and artery, respectively. This model assumes 
that blood enters tissue at the arterial temperature and leaves 
at an equilibrated temperature which is the same as the local 
tissue temperature, regardless of the tissue vasculature. This 
assumption has been questioned by many studies in the lit- 
erature. It has been shown that thermal equilibration between 
blood and tissue occurs before the blood reaches the capil- 
lary bed, and heat exchange between blood and tissue in 



aorta and large arteries (D > 3 mm) is negligible in terms of 
affecting the blood temperature [294,295]. Further theoreti- 
cal studies have shown that thermal equilibration between 
blood and tissue actually occurs in precapillary arteries (or 
arterioles) that are approximately 50 jum in diameter [296]. 
On the venous side, blood temperature begins to deviate 
from its surrounding tissue in approximately the terminal 
vein class (or venules also ~ 50 um in diameter) [296]. Blood 
vessels with a diameter between 50 jum and 3 mm are there- 
fore termed thermally significant vessels. These results show 
that the heat transfer process between blood and tissue is 
very different from the mass transfer process. The latter oc- 
curs dominantly in the capillary bed. Therefore, it is impor- 
tant to account for the effect on heat transfer of different 
classes (or diameters) of blood vessels differently for accu- 
rate thermal modeling. The model for blood perfusion in the 
classical Pennes bioheat equation seems inadequate for accu- 
rate thermal modeling. 

Further investigation have shown that arterioles and ve- 
nules that lead to and from the capillary bed lie in closely 
matched pairs with a countercurrent flow pattern, and mod- 
els are therefore developed to take into account of this ob- 
servation [296-305]. In this case, up to three separate but 
coupled governing equations are required to model the cou- 
pled heat transfer processes in arteries, veins, and the tissue 
region as reviewed elsewhere [146,306,307]. Although mod- 
els accounting for the effects of vessel size and countercur- 
rent flow are more theoretically rigorous, their application is 
very limited [146]. This is probably because the more 
mechanistic models require knowledge of the complex anat- 
omy of the vasculature in a specific tissue, which unfortu- 
nately, is often unknown [306]. On the other hand, many 
recent studies have verified the applicability of the Pennes 
bioheat equation in a variety of tissues/organs by experimen- 
tal data, especially in tissues/organs (e.g., kidney cortex) 
where small blood vessels (D < 0.3 mm) dominate [308- 
311]. Blood vessels of 0.3 - 3 mm in diameter could have a 
significant local cooling or heating effect to the adjacent 



64 The Open Biomedical Engineering Journal, 2011, Volume 5 



Xiaoming He 



tissue, which is not considered in the Pennes bioheat equa- 
tion and should be accounted for in thermal modeling with 
special attention. 

With much attention being focused on modeling the 
thermal effect of blood perfusion using more realistic vascu- 
lature data, the magnitude of the blood perfusion (i.e., w) has 
been often taken as constant. However, the blood perfusion 
during heating has been observed to be strongly dependent 
on both heating time and temperature [312-314], which has 
been shown to significantly affect the thermal and injury 
modeling [38,315,316]. For example, the change of blood 
perfusion has been modeled as follows [317]: 

w = w 31c e~ n (29) 

where w 37C is the baseline blood perfusion rate at 37 °C that 
can be found in the literature for many types of tissues [146] 
and Q is injury to the vasculature which can be calculated 
using Equation 6. Therefore, this model basically takes the 
blood perfusion change in response to heat as a first order 
kinetic process. This model may be sufficient to predict the 
change of blood perfusion in tumor, which has been ob- 
served to decrease monotonically with the increase of ther- 
mal dose (i.e., the accumulation of heat with time reflected 
by injury, Q) [312-314]. However, it is inadequate to model 
blood perfusion response to heating in normal tissue, which 
has been observed to increase initially followed by a gradual 
decrease when the thermal dose is accumulated to beyond a 
threshold level to result in stasis in the tissue vasculature 
[38]. The initial increase of blood perfusion in response to 
heating in normal tissue has been attributed to the host medi- 
ated physiological response, in which elevated blood perfu- 
sion is used to dissipate heat to bring the temperature back to 
normal. Tumor tissue vasculature is believed to be generally 




Degree of vascular stasis 
predicted using 1 st order model 

Fig. (15). Change of normalized blood perfusion (w/w 37 c) in nor- 
mal porcine kidney tissue during heating as a function of vascular 
stasis predicted using the 1 st order kinetic model: Blood perfusion 
increases initially when the predicted vascular stasis is less than ~ 
10% followed by a gradual decrease and eventually complete stasis, 
indicating the first order model is insufficient in predicting the 
change of blood perfusion in response to heating in normal tissue. It 
is reported that the initial increase in blood perfusion from the base- 
line value is absent in heating tumor tissue in vivo. Figure reprinted 
and redrawn from reference [38] with permission from Informa. 



less well organized and leaky and therefore, lacks such a host 
mediated physiological response [318-320]. To model the 
blood perfusion response to heating in normal tissue, Equa- 
tion 29 can be corrected as follows [38]: 

w = w 37C f(e- Q )e- a (30) 

where f(e" Q ) is a correction function. The blood perfusion 
response in thermal therapy of normal porcine kidney tissue 
has been reported (see Fig. 15) and can be fit using Equation 
30 with the following correction function: 

_ f(33.6x 4 +37x 3 -51.8x 2 +32.2x-7.5)/10 4 ,x > 0.9 
X NPKT [276.9x 4 -315.3x 3 +169.2x 2 -42.9x + 5.1,x<0.9 

(31) 

where the subscript NPKT represents normal porcine kidney 
tissue and x represents exp(-Q) with Q given by Equation 6. 
The two kinetic parameters (AE, and A) in Equation 6 for 
calculating x are 98.73 kJ/mole and 1.81 xlO 1 s" 1 , respec- 
tively [38]. Further studies to quantify the two kinetic pa- 
rameters and determine the correction function (for normal 
tissues only) for various types of living systems are neces- 
sary for the wide application of this time and temperature 
dependent model for blood perfusion and therefore, accurate 
modeling of the thermal response in heating tumor and the 
surrounding normal living tissue. 

4.4. Latent Heat and Its Absorption/Release Pattern 

A major component of biological materials is water, 
which undergoes phase change from liquid to either vapor 
when boiling occurs as a result of heating or ice when solidi- 
fication/freezing occurs as a result of cooling. For an ex- 
tremely slow cooling/heating process, phase change of pure 
liquid water to either the vapor or ice phase should occur at 
its equilibrium freezing (~ 0 °C) or boiling (~ 100 °C) point 
accompanying with a significant amount of release (for 
freezing) or absorption (for boiling) of latent heat, respec- 
tively. Due to the existence of impurities such as salts and 
cryoprotectants in tissue water, the latent heat release during 
freezing of biological materials, however, occurs over a tem- 
perature range called the mushy zone [321]. For thermal 
modeling, the upper boundary (0 ub ) of the mushy zone are 
typically taken as the equilibrium freezing point (9 e f P ) of a 
solution, which is a colligative quantity and can be estimated 
based on the osmolality of the solution as follows [322]: 

0 efp =-l.%5Wsm (32) 

where Osm represents osmolality and the subscript efp rep- 
resents equilibrium freezing point. Osmolality is the ratio of 
the total amount (in mole) of particles (ions and non- 
dissociated molecules) in a solution to the total amount (in 
kg) of solvent (water) in the solution. For a dilute, aqueous 
solution at room temperature (when the density of water is ~ 
1 kg/L), the value of osmolality is similar to that of osmolar- 
ity that is the ratio of the total amount (in mole) of particles 
in a solution to the total volume (in liter or L) of the solution. 
For a typical physiological solution with an osmolality of 
0.3, the equilibrium freezing point is around -0.56 °C accord- 
ing to the equation. The lower boundary (Gib) of the mushy 
zone used for thermal modeling varies significantly from 
studies to studies (typically between -4 and -10 °C and can 
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be as low as -40 °C). Various latent heat release patterns (a) 
in the mushy zone have been used as summarized in [142]. A 
simple but quite frequently used linear release pattern of 
latent heat with regard to temperature is given as follows: 

a<0) = (e - e ub )i(e lb - e ub ), e lb < e < e ub (33) 

It has been reported that the magnitude of latent heat (L) 
has a much more significant impact on thermal modeling 
than its release pattern [323]. Biological samples also con- 
tain lipids, proteins, and other organic and inorganic sub- 
stances with essentially negligible latent heat release com- 
paring with water during freezing. Moreover, some amount 
of water may bind with solutes and never become frozen 
(i.e., the so-called bound water) during cooling. Therefore, 
the latent heat of a biological sample should be less than that 
of pure water that is confirmed by recent studies on measur- 
ing the latent heat of some biological solutions using DSC 
[324,325]. In the absence of experimental data, however, the 
latent heat of a biological sample may be estimated as a first 
cut by using the total tissue water content and the latent heat 
of pure water as follows [142]: 

L = <t>L w (34) 

where (p is total water content in a biological sample in 
weight and the latent heat of solidification for pure water 
(L w ) is 335 J/g. Some studies approximates the latent heat of 
biological tissue as the product of the total water content ((p) 
and the latent heat of a physiological solution (lx PBS, ~ 
300 J/g [324,325]), which is much closer to tissue than pure 
water in terms of chemical composition. The total water con- 
tent of a number of tissues can be found in the literature 
[146,326] and/or determined quantitatively by drying (i.e., 
dry weight vs. wet weight) or magnetic resonance spec- 
trometry. 

Very few studies have been reported on the absorption 
pattern of latent heat during boiling in heating biological 
materials, probably because boiling is typically avoided in 
many thermal therapy protocols. The rationale is that boiling 
off tissue water at a temperature higher than 100 °C leads to 
tissue carbonization, a type of wound or lesion that is some- 
times very difficult to heal in the host. In the event that 
evaporation does occur, the amount of latent heat release 
may be calculated similarly to Equation 34 based on the 
weight fraction of water in the sample. The latent heat of 
evaporation for pure water is 2,250 J/g. Similarly, the pattern 
of latent heat absorption may be modeled as being linearly 
dependent on temperature over a small temperature range 
(e.g., 99- 100 °C) [316]. 

4.5. Thermophysical Properties 

The knowledge of thermophysical properties (mainly 
thermal conductivity and specific heat) of biological tissues 
is crucial for accurate modeling of the thermal response in a 
bioheat transfer process. Many studies have been performed 
to measure the properties of various biological tissues, and 
the results from these studies have been reviewed and tabu- 
lated in the literature [146,285,326,327]. Nonetheless, data 
on thermophysical properties of tissue below -30 °C or 
above 40 °C are still scant. Such data are undoubtedly in 
demand because these temperatures are typically encoun- 
tered in thermal destruction and stabilization applications. 



The main factors that affect tissue thermophysical prop- 
erties are temperature and sample composition (mainly water 
content). Because thermophysical properties of water (> 0 
°C) and ice (< 0 °C) over a wide temperature range (from 
liquid nitrogen to boiling water) are well-known and the total 
water content of many tissues are available 
[38,228,321,326,328], the properties of biological samples 
may be estimated using those of water/ice based on their 
total water content. By lumping together the effects of all 
non-aqueous contents, the specific heat (c) of unfrozen bio- 
materials may be estimated as follows [142]: 

c„ =[0 + 0.4(1 -0)]c w (35) 

where the subscript uf represents unfrozen and the first and 
second terms in the brackets represent the effect on specific 
heat of aqueous and non-aqueous contents in tissue, respec- 
tively. A similar relationship between water content and 
thermal conductivity (k T ) of unfrozen biomaterials is as fol- 
lows [142]: 

^ = [0 + 0.28(1-0)]*^ (36) 

where the first and second terms in the brackets represent the 
effect on thermal conductivity of the aqueous and non- 
aqueous contents in tissue, respectively. One can suggest a 
similar approach based on water content to estimate the 
thermophysical properties of frozen tissue. In this case, 
Equations 35 and 36 can be modified to predict the specific 
heat and thermal conductivity of frozen tissue as follows 
[142]: 

^=0^ + 0.4(1-0)^ (37) 

*7\/=#r,i«+ 0.28(1 -0)*^ (38) 

where the subscript f represents frozen and the first and sec- 
ond terms in Equations 37 and 38 represent the effect of tis- 
sue water and non-aqueous contents on the thermophysical 
properties, respectively. Because the non-aqueous contents 
do not undergo phase change during freezing, their contribu- 
tions to the specific heat and thermal conductivity repre- 
sented by the second term of Equations 37 and 38 are calcu- 
lated based on that of water rather than ice. In addition, one 
might use the thermal conductivity and specific heat of water 
at 0 °C to approximate that of sub- cooled water below 0 °C 
as the first cut considering that the latter are still not well- 
established. Finally, determining the thermophysical proper- 
ties of subcooled and vitrified water or solutions below their 
freezing point is crucial for accurate thermal modeling of the 
bioheat transfer processes in thermal stabilization by vitrifi- 
cation (without ice formation). 

5. SUMMARY AND OUTLOOK 

Temperatures over a wide range (from cryogenic to hy- 
perthermic regimens) have been utilized to either selectively 
destroy unwanted or stabilize important biologicals for the 
treatment of diseases. Thermal destruction can be realized by 
freezing biologicals to below -20 °C using cryoprobes or 
heating them to more than 50 °C using RF/microwave, laser, 
and HIFU as the energy source in a minimally invasive man- 
ner. Cryopreservation of biologicals by slow-freezing is as- 
sociated with injury due to ice formation and freeze concen- 
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tration induced cell dehydration. Conventional vitrification is 
associated with cytotoxicity of high concentration of cryo- 
protectants. Research focus in cryopreservation is shifting to 
achieve vitrification using a low concentration of cryoprotec- 
tants (low-CPA vitrification), which combines the benefits of 
both slow-freezing and conventional vitrification while 
avoiding their shortcomings. Eventually, thermal stabiliza- 
tion at ambient temperature (dry or lyopreservation) is the 
solution to the increasing demand of wide and convenient 
distribution of biopreserved products for modern cell-based 
medicine. However, lyopreservation of eukaryotic mammal- 
ian cells at ambient temperature is still challenging today. 
Integrating modern micro and nanotechnology into the field 
of cell preservation might be the solution to the challenges. 
Studies have been performed to understand both hyperther- 
mic and cryothermic injury at the molecular, cell, and tissue 
levels. The mechanism of injury to biologicals of hyperther- 
mic temperatures is better understood than that of cryogenic 
temperatures. Interestingly, using therapeutic agents such as 
TNF-oc and anticancer drugs to augment thermal injury and 
the targeted/controlled delivery of the agents to tumor site 
using nanotechnology are attracting more and more atten- 
tion. Mechanistic models for hyperthermic injury have been 
developed with success whereas developing such models for 
cryothermic injury is more challenging and very much in 
need. Much more work is still in need for thermal modeling 
of bioheat transfer. For hyperthermic applications, it is par- 
ticularly challenging to model the change of blood perfusion 
with regard to both temperature and time (or the extent of 
vascular stasis) for various types of tissues. For cryogenic 
applications, the heat transfer process is affected by addi- 
tional thermophysical events such as ice formation and glass 
transition. Therefore, measurement of latent heat (as a result 
of ice formation and glass transition) and thermophysical 
properties of various biological materials at cryogenic (either 
frozen or vitrified) temperatures is of great importance for 
accurate thermal modeling of the bioheat transfer process. 
With the development of advanced models strengthened by 
further experimental investigations, the thermostability of 
biologicals over a wide temperature range (from that of liq- 
uid nitrogen to that of boiling water) can be engineered for 
either effective destruction of unwanted biomaterials (e.g., 
tumor) or banking of desired biologicals (e.g., stem cells) for 
the emerging cell-based medicine. 
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